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We argue on the basis of empirical data that Kolmogorov's refined similarity hypothesis (RSH) 
needs to be modified for transverse velocity increments, and propose an alternative. In this new 
form, transverse velocity increments bear the same relation to locally averaged enstrophy (squared 
vorticity) as longitudinal velocity increments bear in RSH to locally averaged dissipation. We sup- 
port this hypothesis by analyzing high-resolution numerical simulation data for isotropic turbulence. 
RSH and its proposed modification for transverse velocity increments (RSHT) appear to represent 
two independent scaling groups. 
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In all small-scale turbulence research driven by expec- 
tations of universality, one deals with so-called velocity 
increments, which are differences of velocities between 
two spatial positions separated by a fixed distance. Of 
special interest are velocity increments when the separa- 
tion distances belong to the inertial range, that is, length 
scales which are small compared to a typical large-scale 
motion but large compared to a typical viscous cut-off 
scale. It is fair to say that much of the phenomenological 
work on this subject, summarized in Refs. is based 

to some degree or another on two sets of similarity hy- 
otheses proposed by Kolmogorov, abbreviated as K41 
and K62 ||. The latter, also called the refined simi- 
larity hypotheses (RSH), have been verified both in real 
experiments and numerical simulations of turbulence — 
at least to an extent that certifies them as reasonable. 
The verification has focused, largely for historical rea- 
sons of experimental convenience, on longitudinal veloc- 
ity increments, that is, velocity increments for which the 
separation distance is aligned with the velocity compo- 
nent considered. The theme of this Letter is that RSH is 
inadequate for transverse velocity increments, for which 
the separation distance is transverse to the direction of 
the velocity component considered, and that a non-trivial 
modification (RSHT) is necessary to account for exper- 
imental facts. We conclude that RSH and RSHT form 
two independent scaling groups in small-scale turbulence. 

We need a few definitions before proceeding fur- 
ther. The rate of energy dissipation is given by e — 
^(dui/dxj + duj/dxi) 2 , where v is the fluid viscosity, 
and its local average is defined as e r (x, t) = j v e dV, 
where the integration volume V(r), with the character- 
istic linear dimension r, is centered at the spatial posi- 
tion x. The longitudinal velocity increment is defined as 
Su r = u(x + r) — u(x), where the velocity component u 
and the separation distance r are both in the same di- 
rection, say x. The form of RSH that has been verified 
extensively p-S is given by the relation 



8u r =j3 1 (re r ) 1 / 3 , 



(1) 



where [3\ is a stochastic variable independent of r and 
e r . Assume now that the p-th order longitudinal struc- 
ture function Sp{r) = ((Su r ) p ) ~ and ~ r Tp . 
Here the angular brackets denote suitable averages. It 
then follows from Eq. (|l|) that 
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(2) 



This equation connects the scaling exponents of the longi- 
tudinal structure functions with those of the locally aver- 
aged dissipation function. If a further assumption about 
the statistics of e r can be made, such as log- normal [0, 
multifractal j|, or log-Poisson po[| , r p can be obtained 
analytically via Eq. (j^) . 

The thinking in K62 is that Eq. (1) holds equally well if 
the longitudinal velocity increment is replaced by a trans- 
verse velocity increment, namely Sv r = v(x + r) — v(x), 
where the separation distance r is transverse to the veloc- 
ity component v. If true, this would imply that the scal- 
ing exponents £j and , for longitudinal and transverse 
structure functions, respectively, are equal. A kinematic 
constraint from isotropy 0,11] assures us that = 
Two sets of measurements [12| appear to suggest (or im- 
ply) that the equality holds, within experimental uncer- 
tainty, for larger p as well. On the other hand, more 
recent numerical [ pL 3|Jl 4f| and experimental work |lq | has 
revealed that transverse velocity increments are more in- 
termittent than longitudinal ones, and that are mea- 
surably smaller at least for p > 4. These results 
imply that the RSH cannot be equally true for both lon- 
gitudinal and transverse velocity increments. With this 
in mind, we first obtain the scaling exponents of longi- 
tudinal and transverse structure functions. We employ 
the numerical data from simulations of isotropic turbu- 
lence, carried out using 512 3 mesh points in a periodic 
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box; a statistically steady state was obtained by forcing 
low Fourier modes. For details see the paper cited first in 
||. The Taylor microscale Reynolds number R\ = 216, 
which is close to the limit of the current computational 
capability. Averages over 10 large-eddy turnover times 
were used in the data analysis. 

Figure 1 shows the flatnesses of longitudinal and trans- 
verse velocity increments as functions of r. For all values 
of r except perhaps those comparable to the box size, 
the velocity increment along the transverse direction has 
a larger flatness, suggesting that the transverse velocity 
increment is more intermittent 14 1. This shows that sig- 
nificantly more averaging is required for transverse quan- 
tities than for longitudinal quantities. 
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FIG. 1. Flatnesses of velocity increments, S% (r) / S2 (r) 2 and 
Sj (r) I S2 \r) 2 , as functions of r. 

In Fig. 2, the transverse structure functions S'J(r) are 
plotted as functions of r for p = 2, 4, 6, 8, 10. An inertial- 
range power-law scaling can be identified for 0.2 < r < 
0.6 (the whole box size being 2ir). This is also the scaling 
range determined [|l6| from Kolmogorov's 4/5 law for the 
third-order structure function |17|. 
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FIG. 2. The transverse structure functions Sp(r) as functions 
of r for p = 2, 4, • • •, 10. 

The exponents determined by least-square fits 
within the inertial range just mentioned, are plotted 



in Fig. 3. The use of Extended Self-Similarity (ESS) 
technique Jl8| yields very similar results |l9| . The re- 
sults for Cp , obtained earlier pc|[ , have also been plot- 
ted for comparison. These results demonstrate that 
while both exponents show considerable anomaly due 
to intermittency effects (as evidenced by the deviation 
from the dotted line given by K41), the transverse ex- 
ponents are systematically smaller than the longitudinal 
ones for p > 3. Typical scaling exponents for Sj are: 
Cj = 0.71 ± 0.04, (I = 1-25 ± 0.067, (J = 1-63 ± 0.079 
and — 1-87 ±0.078. While the longitudinal exponents 
agree quite well with existing scaling models J9], ^0| , ^l| , 
the difference between the transverse exponents and the 
models increases with p. For brevity, we show only a com- 
parison with the log-Poisson model IIM, which is typical. 




FIG. 3. Numerical results for the transverse scaling expo- 
nents, ("J and , as functions of p. The dotted line is for the 
normal scaling relation (K41) and the dashed line is for the 
log-Poisson model fll~c| . 

If the scaling exponents for longitudinal and transverse 
structure functions are different, it is clear that Eq. (1), 
now known to be valid for the former, cannot be valid for 
the latter; some modification will therefore be required. 
One possibility suggests itself from kinematic considera- 
tions. Recall from the second reference of M that RSH 
may be expected because the longitudinal velocity incre- 
ment 5u r is an integral over the separation distance r 
of the velocity derivative du/dx, whose square is one of 
the components of the energy dissipation. In a similar 
way, the transverse structure function Sv r is the integral 
of dv/dx, which is one of the two terms of the vortic- 
ity component oj z = dv/dx — du/dy. It is reasonable to 
guess, then, that a suitable modification of Eq. (1) would 
be 

Sv r = /? 2 (rttr) 1/3 , (3) 

where ft — voj 2 , u> — V x u is the vorticity, and f2 r is the 
local average of in the same way as e r is the local aver- 
age of e in Eq. (Q). A necessary condition for the above 
equation to be plausible is that e r should scale differently 
from £l r and that the difference must be compatible in 
sign with that of the difference between and £j . That 
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this is indeed so was demonstrated recently H); that is, 
if (QP) ~ r°p, o p < T p . Thus, Eq. (3) is a priori worth 
exploring as a candidate for the refined similarity hypoth- 
esis for transverse velocity increments (RSHT) . 

We now test the validity of Eq. (3). In particular, we 
test whether the relation 



t T = P 



Op/3, 



(4) 



which follows from Eq. (3), is true. We should note that 
the equations 



Sv r =fo{re r ) 1 l\ 



(5) 



Su r = (3 4 {rfl r ) 1/3 (6) 
are both dimensionally plausible, and would yield 



p 

3 



Op/3, 
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(8) 



Equation (4) would derive greater support if the alterna- 
tive relations (7) and (8) can be shown to be unsatisfac- 
tory. 




FIG. 4. Verification of the scaling relations in RSHT. (a) 
shows the results for the transverse structure functions and 
(b) for the longitudinal case. The dashed lines correspond to 
p/3, which would be followed by the numerical data if the 
similarity hypothesis used is correct. See the equation pairs 
(2) and (4), and (7) and (8). 

In Figs. 4 (a) and (b), we show comparisons of numeri- 
cal results for p up to 10 with predictions of (|j) and (||) for 
the transverse case, and with predictions of (Q) and (0) 
for the longitudinal case. (For odd order, we have used 



absolute values of the appropriate velocity differences.) 
The exponents r p / 3 and o p /3 used in the above relations 
were taken from previously published results for the same 
data [^o|,^2). For p < 3, the differences between the im- 
plied relations (|2|), (||), and (Q), (||) and numerical results 
are relatively small, and cannot be distinguished easily. 
On the other hand, with increasing p, only £ p — Tp/3 and 
(T — Op/3 coincide with p/3, showing that Eqs. (||) and 
(||) are valid to a good approximation. The numerical 
evidence that £p — o p / 3 and Cj — T P /3 depart from p/3 
for large p effectively disqualifies (||) and (|J) . Note that 
the error bars in Figs. 4 are quite small and do not affect 
this conclusion. 

A more detailed test of Eq. (^) is given in Fig. 5, 
where we plot the normalized probability density func- 
tion (PDF) of /3 2 = Av r /(rfl r y/ 3 with r = 0.2, 0.34 
and 0.44, all of which lie in the inertial range. The three 
PDFs collapse quite well. It is seen that in the region 
(7) /V(/3f) 1/2 <3, the PDFs agree with the standard Gaus- 
sian distribution. For larger values of #2/ ((/?2) 2 ) 1 ^ 2 , the 
PDFs seem to depart slightly from Gaussian, with a ten- 
dency perhaps towards exponential. This trend has also 
been seen for the longitudinal counterpart, investi- 
gated by Wang et al. ||. One difference between the 
PDFs of Pi and fa, which follows from Kohnogorov's 
4/5-ths law Jl7]], is that /?2 is not expected to be skewed, 
whereas (3\ should have a third moment equal to —4/5. 
The PDFs of {3$ and /?4 have also been obtained but not 
shown here. They do not collapse for inertial range sepa- 
rations; nor do they agree with the Gaussian distribution 
even in the core region; and they possess larger depar- 
tures from Gaussian near the tails. 

The results in Figs. 4 and Fig. 5 suggest that RSHT 
proposed in (|^) is a good working approximation con- 
necting the transverse velocity increment with the en- 
strophy field. In addition, the differences of the intermit- 
tency properties between 5u r and Sv r as well as e r and 
fl r strongly suggest the possibility of two independent 
scaling groups. Note that e r and fl r are the symmetric 
and antisymmetric parts of the strain-rate tensor, respec- 
tively; so we expect that this difference, reflecting the 
difference in the inertial-range physics of the dissipation 
and vortex dynamics, may quite plausibly carry over to 
the high Reynolds number limit as well. 

It should be pointed out that a more general refined 
similarity hypothesis, which encompasses results for both 
velocity components, cannot be ruled out. One possible 
scenario could be 



Av r = ftWeJ-") 1 / 3 , 



(9) 



where /3g is again a universal stochastic variable and a 
could, in general, be a function of the order index p [ p3[ . 
Equation ^ is more complex than ([!]) and (||) , since the 
scalings of the cross-correlation functions (^ e r) have to 
be taken into account. 
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FIG. 5. Normalized PDFs of /3 2 /((/3 2 ) 2 > 1/2 for r = 0.2, 0.34 
and 0.44. 

In summary, we have proposed a refined similarity 
hypothesis for the transverse direction (RSHT); it con- 
nects the statistics of the transverse structure function 
in the inertial range with the locally averaged enstrophy. 
Data analysis of DNS for the Navier-Stokes equations at 
moderate Reynolds numbers demonstrate that RSHT is 
valid for the transverse structure functions in the inertial 
range, while Kolmogorov's RSH has been shown previ- 
ously to be essentially correct for the longitudinal case 
p-p|; the latter result is confirmed again in this paper. 
The important implication of RSHT is the possibility of 
the existence of two independent scaling groups which 
may correspond to different intermittent physics in fluid 
turbulence: these are related to the symmetric part of 
the strain-rate, or the dissipation physics, and those re- 
lated to the antisymmetric part of the strain-rate, or the 
vortex dynamics. This view also suggests that no more 
than two sets of independent exponents are required for 
describing the scaling of all small-scale features as longi- 
tudinal and transverse velocity increments, dissipation, 
enstrophy, and circulation M]. 
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